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Abstract — We study a four terminal parallel relay- 
eavesdropper channel which consists of multiple independent 
relay-eavesdropper channels as subchannels. For the discrete 
memoryless case, we establish inner and outer bounds on the 
rate-equivocation region. For each subchannel, secure trans- 
mission is obtained through one of the two coding schemes 
at the relay: decoding-and-forwarding the source message or 
confusing the eavesdropper through noise injection. The inner 
bound allows relay mode selection. For the Gaussian model 
we establish lower and upper bounds on the perfect secrecy 
rate. We show that the bounds meet in some special cases, 
including when the relay does not hear the source. We illustrate 
the analytical results through some numerical examples. 

I. Introduction 

The relay channel has been analyzed in (TJ, J2 (and 
references therein), but the focus was on how to increase 
achievable rate and reliability. The idea of cooperation be- 
tween users in the context of security was introduced in J5] 
(and references therein). The premise is that when the main 
channel is more noisy than the channel to the eavesdropper, 
cooperation between users is utilized to obtain a positive 
secrecy capacity. Secrecy is achieved by using the relay as 
a trusted node which facilitates the information decoding at 
the destination while confusing the eavesdropper. 

In this paper, we study a parallel relay-eavesdropper chan- 
nel. A parallel relay-eavesdropper channel consists of a gen- 
eralization of the setup in [3 1 to the case in which each of the 
source-to-relay (S-R), source-to-destination (S-D), source-to- 
eavesdropper (S-E), relay-to-destination (R-D) and relay-to- 
eavesdropper (R-E) link is composed of several independent 
parallel channels as subchannels. The model is depicted in 
Figure Q] For this model, we establish outer and inner bounds 
on the rate-equivocation region for the discrete memoryless 
case. The inner bound is obtained with a coding scheme in 
which, for each subchannel, the relay operates in decode- 
and-forward (DF) or noise forwarding (NF) mode. The outer 
bound does not follow directly from the single-letter outer 
bound for the relay-eavesdropper channel developed in J3] 
Theorem 1] and so, a converse is needed. This converse 
includes a redefinition of the involved auxiliary random 
variables, a technique much similar to the one used before in 
the context of secure transmission over broadcast channels 
E]. We also show that the bounds on the equivocation 
rate coincide in the case in which all the subchannels are 
degraded, thus characterizing the secrecy capacity. 

For the Gaussian model, we focus on the perfect secrecy 
case. We establish lower and upper bounds on the perfect 
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Fig. 1. Parallel relay-eavesdropper channel. 



secrecy rate. We note that establishing a computable upper 
bounds on the secrecy rate for the Gaussian model is non- 
trivial. In part, this is because converse techniques that 
are obtained directly from the analysis in the discrete case 
involves auxiliary random variables, the optimal choice of 
which is difficult to obtain. We develop a new upper bound on 
the secrecy rate for the parallel Gaussian relay-eavesdropper 
channel. Our converse proof uses elements from converse 
techniques developed in Q, J6) in the context of multi- 
antennas wiretap channel; and in a sense, can be viewed 
as a partial extension of these results to the case of the 
studied model. The established upper bound on the secrecy 
rate shows some degree of separability for different parallel 
subchannels. It is especially useful when the multiple access 
part of the channel is the bottleneck. 

We also study a special Gaussian case in which the 
relay does not hear to the source, for example due to very 
noisy source-to-relay links. In this case, we show that noise- 
forwarding on all links achieves the secrecy capacity. The 
converse proof follows from the general converse established 
for the general Gaussian case, and a new genie-aided upper 
bound that assumes full cooperation between the relay and 
the destination, and a constrained eavesdropper. The eaves- 
dropper is constrained in the sense that it has to treat the 
relay's transmission as unknown noise for all subchannels, 



an idea used previously in the context of a class of classic 
relay-eavesdropper channel with orthogonal components in 
J7J. These assumptions turn the parallel Gaussian relay- 
eavesdropper channel into a parallel Gaussian wiretap chan- 
nel, the secrecy capacity of which is established in J4). 

Furthermore, we apply the results developed for the paral- 
lel Gaussian relay-eavesdropper channel to the fading relay- 
eavesdropper channel, which is a special case of parallel 
Gaussian relay-eavesdropper channel, with each fading state 
corresponding to one subchannel. We illustrate our results 
through some numerical examples. 

In this paper, the notation -X'ml] is used as a shorthand 
for (X\,X 2 , ■ . ■ , X£), the notation X"^ is used as a 
shorthand for (Xf , X£, • • • , X£) where for I = 1, . . . , L, 
X\ l := (Xn,Xi2, ■ ■ ■ ,Xi n ) and the notation <^i[i,l] is use d 
as a shorthand for Xu X X 12 . . . x X\l. We define the function 
C(x) as | log 2 (l + x). Throughout the paper the logarithm 
function is taken to the base 2. 

II. The parallel Relay-Eavesdropper Channel 
A parallel relay-eavesdropper channel is a four terminal 
network consisting of X^ix] , ^2[i,L] as finite input alphabets 
and y[L l L],yi[i,L]>y2[x,L] as finite output alphabets. The 
transition probability distribution is given by 



J L P(yi,yu,V2l I xu,x 2 i) 



(1) 



i=i 



where x u £ X iu x 2l £ X 2U y u £ yu,yi £ y t and y 2l £ y 2 i, 
for I = 1, • • • , L. For the subchannel I, Xu and X 2 i are the 
inputs from the source and relay; and Yu , Y 2 i , Yi are the out- 
puts at the relay, eavesdropper and destination, respectively. 
The source sends a message W € W = {1, • • • , 2 nR } using 
a (2 nR ,n) code consisting of: 1) a stochastic encoder at the 
source that maps W — s- X™, t L y 2) a relay encoder that maps 

/*(^i[i,L]) ~^ X 2[i,L],i for 1 < i < n, and 3) a decoding 
function giY[i v) ~^ W. The average error probability of a 
(2 nR , n) code is defined as 



P e =iE P{9(Y{l L] ) ^ W\W}. 



(2) 
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The eavesdropper listens to what the source and relay trans- 
mit for free, due to the wireless nature of the medium. It then 
tries to guess the information being transmitted. Denoting 
^2?i l] me 0Ut P ut at tne eavesdropper, the equivocation rate 
per channel use is defined as R e = H(W\Y£ 1 ix)/n. Perfect 
secrecy for the channel is obtained when the eavesdropper 
gets no information about W from Y^L L , . That is, the 
equivocation rate is equal to the unconditional source entropy. 
A rate equivocation pair (R, R e ) is achievable, if for any 
e > there exists a sequence of codes (2 nR , n) such that for 
any n > n(e) 
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III. Discrete memoryless channel 

In this section, we establish outer and inner bounds on 
the rate-equivocation region of a parallel relay-eavesdropper 
channel. 

A. Outer bound 

Theorem 1: For the parallel relay-eavesdropper channel 
with L subchannels, and for any achievable rate-equivocation 
pair (R, R e ), there exists a set of random variables Ui 

(Vu,V 2l ) -> (X lh X 2l ) -> (Y h Y lh Y 2l ), l = l,...,L, such 
that (R, R e ) satisfies 

R < win(^2l(V 1 iV 2 i;Y l ),J2l(V ll ;Y l Y ll \ V 2l )\ 
^ i=i i=i > 

R e < R 

Re < miii [ KYuVmY, I U t ) - I(VuV 2l ;Y 2l \U,), 
^ i=i 

^HVuiYYu I V 2 iUi) - I(VuV 2l ;Y 2l \ U t )\. (4) 
i=i ' 

Proof: The proof of Theorem 1 is given in Appendix Q] 

□ 

B. Achievable rate-equivocation region 

Theorem 2: For the parallel relay-eavesdropper channel 
with L subchannels, the rate pairs in the closure of the convex 
hull of all (R, R e ) satisfying 



R < mini ^IiYuV^Y^UO^IiVuiYulViiUi) 
^ leA 

+ I{Vu;Yi\V 2 i) 



leA 



l£A c 



Re < R 



Re < mm{Y J I { v nV2V,Y l \Ui)-I{ViiV 2 r,Y 2l \Ui), 
^ leA 

Y HVu ; Y u | VuUi) - 1( VuV 2 i; Y 2l \ U x ) X 
leA > 

+ Y I(Yu',Yi\V2i)+nan( Y ItyaM), 
ieA c ^ ieA c 

Y I(V2i;Y 2l \Vu)\-mmij2 T ^ Y ^ 
ieA c J ^ ieA c 

Y l(V2i;Y 2l )\ - Y l(Vu;Y 2l \V 2t 



(3) 



^21), (5) 
ieA c ' ieA c 

for some distribution p(ui,vu,v 2 i,xu,x 2 i,yi,yu,y 2 i) = 
p(ui)p{vu,v 2l \ui)p(xu,x 2 i\vii,v 2 i)p(yi,yxi,y 2 i\x 1 i,x 2 i) 
for I £ A and p(vu,v 2 i,xu,x 2 i,yi,yii,y 2 i) = 
p{v u )p{v 2 i)p(xii\vii)p{x 2 i\v 2 i)p(yi 1 y 1 i 1 y 2 i\x 1 i,x 2 i) for 
I £ A c , are achievable. 

In the statement of Theorem 2, sets A and A c represent the 
subchannels for which relay operates in DF and NF mode, 
respectively. The region in Theorem |2] is obtained through a 
coding scheme which combines appropriately DF and NF 



schemes. The rates for the DF scheme can be obtained 
readily by setting U := U^^Vi := V 1[hlAl] ,V 2 := 

V 2 [1,\A\]jY : = ^[1,1^1]' Y l : = and Y 2 : = *2[1,|.A|]> 

for I G A in [3] Theorem 2]. Similarly the rates for 
NF scheme can be readily obtained by setting V\ := 



Vl[l,\A<=\],V2 

and Y 2 := i2[i,|X c |]> f° r ' "= -A c in |3, Theorem 3]. 



^[1,1^1], y : = y[i,|^|], y 



y 



l[l,|-4<= 



Remark 1: For a parallel relay-eavesdropper channel in 
which all subchannels are degradecQ, i.e., 

V{VUV\UV2l I xu,x 2 i) 

= P(yil I Xu,X 2 l)p{yi I yihX2l)p{V2l I vi,yu,xu > x2i), 
for I = 1, . . . , L, the perfect secrecy capacity is given by 

C s =maxmin [ ^[/(^y a; y | £/,) - /(^ijVa,; Faj I E/j)]~ 
^ 1=1 



Y^[I{V U ;Y U I v 2i [/ z ) - I{V u V 2l] Y 2l 1 1/,; 



/=1 



(6) 



where maximum is over Ui —> (Vu,V2i) —> (Xu,X 2 i) — s- 
{Y u Y u ,Y 2l ), forl = l,...,L. 

Proof: The achievability follows from Theorem 2 by 
setting A c = 0. The converse follows along the lines of 
Theorem 1 and is omitted for brevity. □ 

IV. Gaussian channel 

In this section we study a parallel Gaussian relay- 
eavesdropper channel. We focus on perfectly secure achiev- 
able rates, i.e., (R,R e )= (R,R). The received signals at the 
relay, destination and eavesdropper are given by 

Yil,i = Xn.i + Zn.i 
Yi,i = Xu^i + \ZpTiX 2 i,i + Zij 
Y 2 Lt = Xu ti + «fp2iX 2 i,i + Z 2 i,i (7) 

where i is the time index, {Zu^}, {Zi^} and {Z 2 i_i} are 
noise processes, independent and identically distributed (i.i.d) 
with the components being zero mean Gaussian random 
variables with variances af,, of and o-\\\ Xu,i and X2i,i 
are the inputs from the source and relay nodes respectively. 
The parameter pu indicates the ratio of the R-D link signal- 
to-noise (SNR) to the S-D link SNR and p 2 i indicates the 
ratio of the R-E link SNR to the S-E link SNR for the I th 
subchannel. The source and relay input sequences are subject 
to the following average transmit power constraints 



L n 



1 " 

-zZzZ E iK,i]<Pu 

1=1 i=l 
1 L n 



(8) 
(9) 



A. Lower bound on the perfect secrecy rate 

For the parallel Gaussian relay-eavesdropper channel, de- 
fined by (01, we apply Theorem [2] to obtain a lower bound 
on the perfect secrecy rate. 

Corollary 1: For the parallel Gaussian relay-eavesdropper 
channel (O, a lower bound on the perfect secrecy rate is 
given by ( fTot . 

Proof: The achievability follows by applying Theorem 
[2] with t he cho ice Ui = constant, Vu=Xu, V2i=X2i, Xu = 
X u + ^Su^LX2i, ai = l- a u Xu ~ Af(0, otiP u ) indepen- 
dent of X 2t ~ Af(0,P 2 i) for I e A; and X u ~ Af(0,Pu) 
independent of X21 ~ A/"(0, P21) for I <G A c . Straightforward 
algebra which is omitted for brevity gives ( fTOb . □ 

The parameters Pu and P21 indicate the source and relay 
'power allocated for transmission over the subchannel I. In 
([Tol l, after some straightforward algebra, the contribution 
to the equivocation of information sent through NF (set 
A c ) can be condensed by observing that we only need 
to consider min{^ ie ^c I(X 2 u Y21), J2ieA a I ( X 2i, Yi)} = 
Y^ieA" I(X2i,Yii) in set A c , to get a higher secrecy rate. 
A simplified expression for i?j, ow is given by ( fTTT i. 



B. Upper bound on the perfect secrecy rate 

The following theorem provides an upper bound on the 
secrecy rate of the parallel Gaussian relay-eavesdropper chan- 
nel. 

Theorem 3: For the parallel Gaussian relay-eavesdropper 
channel (|7]), an upper bound on the secrecy rate is given by 



R7< max Y / I(XuX2i:Y l )-I(X ll X 2 i;Y 2 i) (12) 
{K P( e/Cp,}i = i...i ^— ^ 



where the maximization is over [Xi;,^;] ~ 7V(0,Kp ( ) 



with Kp. 



K P , : Ki 



Pu 1P1 VP11P21 

P21 



1 < 



ijji < 1 \, for I = 1, 



1=1 1=1 



.V/ VP11P21 

L, with the covariance matrices 

E{X 1[hL] Xf [hL] }, E[X 2[lii] X^ 1;L] ] satisfying ® and © 
respectively. 

Proof: The result in Theorem 1 established for the DM 
case can be extended to memoryless channels with discrete 
time and continuous alphabets using standard techniques [5] 
Chapter 7]. Taking the first term of the minimization in the 
bound on the equivocation rate, we get 



Re < max^2l(yuVu]Yi \ U t ) - I{V u V 2 v,Y 2l \ Ui) (13) 
1=1 



'in parallel relay-eavesdropper channel if all the subchannels are de- 
graded, the entire relay-eavesdropper channel may not necessarily be de- 
graded. 



where Ui -> {V xl ,V 2l ) -> {X 1U X 21 ) ->■ (Y U Y 1U Y 21 ), for 
I = 1, . , . , L. The rest of the proof uses elements from related 



works in J4) and I0. Continuing from ( TOl l, we obtain 

L 

Re<Y, HVuVu;Yi | Ui) - I(VuV 2l ;Y 2l \ Ut) 
1=1 

(a) L 

< J2 1{VuV 2l -Yi) - I(VuV 2 i;Y 2l ) 
1=1 

L 

< I(VuV 2l ;Y[Y 2l ) - I{VuV 2l ; Y 2l ) 



^jyiXuX^YtYu) - I(XuX 2l ;YiY 2l \ V u V 2 i)\ 



1=1 

L 



1=1 



-[I{X u X 2l ;Y 2l ) - I{X u X 2l ;Y 2l \ V u V 2l )} 

L 

■ Y^IiXuX^Y^) - I(X n X 2l ;Y 2l )} 



i=i 



-[IiXuX^YYu | V u V 2l ) - I(X n X 2l ;Y 2l | V u V 2l ) 

L 

<^yi{X u X 2l -YiY 2 {) - I(X n X 2l ;Y 2l )} 



i=i 

L 



Y^HXiiX^Y \ Y 2l ), 



(14) 



i=i 



where (a) follows by noticing that I(VuV 2 i;Yi \ Ui) — 
I{VuV 2 i;Y 2 i | Ui) is maximized by setting C/;=constant and 
(6) follows from the Markov chain condition (Vu,V 2 i) — > 
(X lh X 2l ) -> {Y h Y lh Y 2l ), for I = 1, . . . , L. 

We now tighten the upper bound (TBI) by using an argument 
previously used in J5), J6) in the context of multi-antennas 
wiretap channels. More specifically, observing that, the origi- 
nal bound ( TOT l depends on p(yi,y 2 i\xu, x 2 i) only through its 
marginals p{yi\xu,x 2 i) and p(y 2 i\xu, x 2 i), the upper bound 
( TBI can be further tightened as 

L 

Re < max Y min I{X u X 2l ;Y{ \ Y 2 \) (15) 
p(x u ,x 2 i) ^ p(,y',,y' 2l \x u ,x 2 i) 

where the joint conditional p(y' ll y' 2l \xu 1 x 2 i) has the 
same marginals as p{yi,y 2 i\x lh x 2 i), i.e., p(y,'|a;iz, x 2l ) = 



p(?//|a;i/,X2i) and p^l^i/,^;) = p{y 2 i\xu, x 2 i). 

It can be easily shown that the bound in (TT3T > is 
maximized when the inputs are jointly Gaussian, i.e., 

[X u ,X 2l ] ~ AA(0,K P( ), K P , G £ P , with /C P) = (k p , : 



4>i VP11P21 P21 



-1 < ipi < 1 >, for 



1, . . . ,L with the covariance matrices EfX^^jX^ L j] and 
E[X 2 [i,l]-X^[i,l]] satisfying © and © respectively 0, J6). 

Next, using the specified Gaussian inputs, it can be shown 
that the evaluation of the upper bound ( TOl l minimized over 
all possible correlations between Y( ,Y 2 \, for / = 1,...,L 
yields 



Re < 



max V I(X u X M ;Yi) - I(X lt X 2l ; Y 2l ). (16) 

{Kp^KpJti...!^ 



□ 



This concludes the proof. 

We now study the case in which the links S-R are very noisy, 
i.e., the relay does not hear the source. 

Theorem 4: For the model (0 in which the relay does not 
hear the source, the secrecy capacity is given by 



C s = min < max C 



1=1 



Pu + PuP 2 i\ _ c ( Pu + 9%P: 



■21 



'21 



1=1 v « 



Pu 



' 2 l 



PIIP-II 



(17) 



where the maximization is over {Pi;, P 2 i}, for / = 1 . . . L, 
such that Y,i=i p ii < p i and Ef=i p 2i < Pi- 
Proof: 

Upper Bound: The bound given by the first term of the 
minimization in ( TT7T > follows from a straightforward applica- 
tion of the result in Theorem 3 — taking independent source 
and relay inputs since the relay does not hear the source 
transmission in this case. 

The bound given by the second term of the minimization 
in (Tl7] i can be established as follows. Our approach borrows 



T-)l0W 

K e — 



max mm 

T.l = lPu<Pl,T,l = l P 2l<P2, 
0<a,<l, fori = 1, ... , LAI 



E^ 



Pu + P11P21 + 2yJa lPll P u P 2 



Pu + P21P21 + 2y/aiP2lPllP 2 , 



r ( P11P21 



-c 



Pu + P21P21 + 2y/a lP 2iPuP 2, 

P21P21 
2 

J 2l 



n ( PllP2l \ r ( P2lP 

£ c te)h£ c (iO- '"» 



T-jlOW 

n e — 



max mm <; 

T.t =1 Pll<Pl,T.Ll P 2l<P2, Vf^A 
0<a,<l, for! - 1, ... , L4 



Pu + PuP-a + 2^/aipuPuP2 



E 



aiPu 



-C 



-C 



Pu + P21P21 + 2^a iP2l P u P 2 , 



Pu + P21P21 + 2^a lP2l P ll P 2: 



{ E 



, , Pll + PllP2l \ _q( Pll+ P21P2I 



^E 



Cl-^)+c(^L\-Cl Pu+p2lP21 



(11) 



elements from an upper bounding technique that is used 
in 0, and can be seen as an extension of it to the case 
of parallel relay-eavesdropper channels. Assume that all the 
links between the relay and the destination are noiseless, and 
the eavesdropper is constrained to treat the relay's signal 
as unknown noise. As mentioned in J7], any upper bound 
for this model with full relay-destination cooperation and 
constrained eavesdropper, also applies for the general model. 

Now, for the model with full relay-destination cooperation 
and constrained eavesdropper, we develop an upper bound on 
the secrecy capacity as follows. In this case, the destination 
can remove the effect of the relay transmission (which is 
independent from the source transmission as the relay does 
not hear the source), and the equivalent channel to the 
destination can be written as 

Y{ ti = Xu,i + Zi,i. (18) 

For the constrained eavesdropper the relay's transmission acts 
as an interference, with the worst case obtained with Gaussian 
X 2 [\_l] P\- The equivalent output at the eavesdropper in this 
case is given by 

Yki = X ^ + yjpnnxiij + Z*,i- d9) 

The rest of the proof follows by simply observing that the 
resulting model (with the worst case relay transmission to 
the eavesdropper and full relay-destination cooperation) is, in 
fact a parallel Gaussian wiretap channel, the secrecy capacity 
of which is established in J4), 

L 

C s < max^7(X 1( ; Y{) - I(X U ; Y 2 \) (20) 
i=i 

where the maximization is over Xu ~ Af(0, Pi/), X21 ~ 
J\f(0,P 2 i) for I = 1...L, with Ya=i P u < p i and 
£f =1J P 2 /<P 2 . 

Finally straightforwad algebra which is omitted for brevity 
shows that the computation of ( f20l > gives the second term of 
the minimization in (T7\ . 

Lower Bound: The achievability follows by computing 
the lower bound in Theorem 2 with the choices |.A C |:=L, 
Vif.=X u , V 2 i:=X 2 i, and X u ~ Af(0,P u ) independent of 
X 2l ~Af(Q,P 2l ). □ 

V. Fading relay-eavesdropper channel 

We consider a fading relay-eavesdropper channel that is 
corrupted by multiplicative fading gain processes in addition 
to additive white Gaussian noise (AWGN) processes. The 
received signals are given by 

Yi t i = h sri Xii + Zii 

Yi = h s d t iXij + h r d,iX 2; i + Zi 
Y% t i = h se ^X\.i + h re ^X 2 .i + Z 2 ^ (21) 

where i is the time index, h s d,i, h r d,i, h se ,i< h re ,i and 
h sr ,i are the fading gain coefficients associated with S-D, 
R-D, S-E, R-E and S-R links, given by complex Gaus- 
sian random variables with zero mean and unit variance 



respectively. The noise processes {Zi.i}, {Zi}, {Z 2 .i} are 
zero mean i.i.d complex Gaussian random variables with 
variances a\, a 2 and a\ respectively. The source and relay 
input sequences are subject to an average power constraint, 

EHiVmI 2 ] < nPi, EtiM^A 2 } < nP 2 . Let 
hi := [h s d,i Kd,i h se .i h re ,i h sr ,i] and we assume that 
perfect channel state information (CSI) is available at all 
nodes, i.e, each node has access to the instantaneous CSI and 
its statistics. For a given fading state realization hi, the fading 
relay-eavesdropper channel is a Gaussian relay-eavesdropper 
channel. Therefore, for a given channel state with L fading 
state realizations, the fading relay-eavesdropper channel can 
be seen as a parallel Gaussian relay-eavesdropper channel 
with L subchannels. The power allocation vectors at the 
source and relay are denoted by Pi(h) and P2{h) respec- 
tively. The ergodic achievable secrecy rate of the fading 
relay-eavesdropper channel (fSTJ, which follows from (fTTI) 
is given by d22l . The upper bound for the fading relay- 
eavesdropper channel follows from ( fT2l and is given by d23l ). 
In the achievable region we proposed a coding scheme which 
is a combination of DF and NF scheme. A pertinent question 
is how to decide which scheme to use on each subchannel 
? To accomplish this we define A:= {h : \h s d\ 2 < \h sr \ 2 } 
contains all the fading state realizations in {h} where the 
S-R link is better than S-D link. The complement of set A 
is A c := {h : \h sd \ 2 > \h sr \ 2 }. 

VI. Numerical Results 

We consider a fading relay-eavesdropper channel with 
L fading states. It is assumed that perfect channel state 
information is available at all nodes. We can consider this 
channel as a Gaussian relay-eavesdropper channel with L 
subchannels. We model channel gain between node i G {s, r} 
and j £ {r, d, e} as distance dependent Rayleigh fading, that 
is, hi j = h\ J , where 7 is the path loss exponent and 
h'i j is a complex Gaussian random variable with zero mean 
and variance one. Each subchannel is corrupted by AWGN 
with zero mean and variance one. The objective function for 
both lower and upper bounds are optimized numerically using 
AMPL with a commercially available solver, for instance 
SNOPT. Furthermore, for each symbol transmission same 
subchannel is used on S-R and R-D links to make the 
optimization tractable. 

To illustrate the system performance, we set the source and 
relay power to 64 Watt each. We consider a network geometry 
in which the source is located at the point (0,0), the relay is 
located at the point (c?,0), the destination is located at the 
point (1,0) and the eavesdropper is located at the point (0,1), 
where d is the distance between the source and the relay. 
In all numerical results we set path loss exponent j:=2. Fig. 
|2] shows the power allocation for a fading channel with 64 
subchannels where the relay is located at (0.5,0), and marker 
'x' denotes NF on a particular subchannel while marker 
'□' denotes DF on a particular subchannel. It can be seen 
from Fig. [2] that, achievable perfect secrecy rate is zero for 
some subchannels. Roughly speaking, this happens when the 
condition \h r d\ 2 > \h re \ 2 is violated. 
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Fig. 2. Achievable perfect secrecy rate of a fading parallel relay- 
eavesdropper channel. 



Fig. 3. Bounds on perfect secrecy rate. 



Mode selection at the relay by only considering the relative 
strength of the S-D and the S-R link is suboptimal because 
the achievable secrecy rate (|22| | also depends on the other 
link gain. We now consider the case in which the relay 
independently selects the scheme which maximize the rate for 
each subchannel. When the relay is close to the source, it uses 
DF scheme on all subchannels. Similarly when the relay is 
close to the destination, use of NF scheme on all subchannels 
offers better rate. The region when the relay is between 
0.5 < d < 1.2, it selects between DF and NF scheme based 
on link gains of S-D, S-R link as mentioned above. In Fig. [3] 
we plot the optimized lower and upper bounds on the secrecy 
rate for fading relay-eavesdropper channel. It can be seen that 
when the relay is close to the source the lower and upper 
bound matches. This follows because of using DF scheme 
on all subchannels. 



Appendix I 

The proof generalizes the results of Theorem 1 in O and 
uses elements from a similar proof in the context of parallel 
BCC in E). 

a) We first bound the equivocation rate as follows. 
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where e„ -} as n 4 oo; (a) follows from Fano's 
inequality; and (6) and (c) follows from lemma 7 in J9)- 

We introduce a random variable T uniformly 
distributed over {1, 2, •••,n} and set, = 
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(Ui,Vu, V2i,X 1 i,X 2 i,Yi,Y 1 i,Y 2 i) satisfies the following 
Markov Chain condition 
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where (/) and (g) follow from the above definition, 
b) We now bound the rate R as follows. 
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where (d) and (e) follow by using the above definition. 
We can also bound the equivocation rate as follows. We 
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Hence, we have 
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where (/i) follows from Fano's inequality; (i) and (j) follows 
from the fact that conditioning reduces entropy. 
We can also bound the rate R as follows 
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where (k) follows from Fano's inequality; and (I) follows 
from the fact that conditioning reduces the entropy. 

Therefore an outer bound on the achievable rate equivoca- 
tion region is given by the following set: 

|J j (R, R e ) that satisfy 435}, 435), ©, <[32) 1 (33) 



where the union is over all probability distributions 

P(«[l,i]j«l[l,i]j«2[l,£]5^1[l,L],iC2[l,i]jl/[l,Lb_yi[l,i],y2[l,i]). 

Finally we note that the terms in d$6]l, 42SK 43(1 . 
and ([32) depend on the probability distribution 

P(U[l,L],Vl[l,L],V2[l,L],Xl[l,L],X2[l,L],V[l,L],yi[l,L],y2[l,L]) 

only through p(ui,vn,v 2 i,xu,x 2 i,yi,yu,y 2 i). Hence, there 
is no loss of optimality to consider only those distributions 
that have the form 
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This completes the proof of Theorem 1 . 
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